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ABSTRACT 

We present a method for the calculation of the effects of external irradiation on 
the geometrical shape of the secondary in a close binary containing a compact star, the 
source of the radiation, and a normal companion star, where we include the possibility 
of shadowing by an accretion disc. The model is based on a simple modification of 
the standard Roche binary potential in which the radiation-pressure force is parame- 
terised using the ratio of the radiation to the gravitational force. We have constructed 
numerical solutions of approximate 3-dimensional irradiated equipotential surfaces to 
demonstrate the main geometrical effects of external radiation pressure. For systems 
in which the inner Lagrangian point is irradiated directly and for sufficiently high 
irradiation fluxes, the critical condition for which the secondary fills its tidal lobe 
changes from an inner to an outer critical configuration, where the critical equipoten- 
tial surface connects to one of the outer rather than the inner Lagrangian point. Such 
a situation may apply to evaporating binary pulsar systems (e.g. PSR 1957-1-20), stars 
orbiting supermassive black holes in AGN and some high-mass X-ray binaries (in- 
cluding Centaurus X-3). For systems containing an accretion disc, which shadows the 
inner Lagrangian point from the external irradiation, the presence of significant radi- 
ation pressure causes a non-axisymmetric deformation of the stellar surface. This has 
particularly important consequences for low-mass X-ray binaries, for which the X-ray 
luminosities can be close to the Eddington limit. We have calculated modified Roche 
potentials to determine the main effects on the optical lightcurves and radial velocity 
curves for typical binary parameters. Compared to previous studies, the inclusion of 
irradiation-pressure effects results in changes in the derived system parameters (e.g. 
component masses, radial velocities) that may be as high as ^ 30 per cent. We con- 
clude that the proper inclusion of irradiation pressure effects is essential for a reliable 
analysis of close binary systems in which the secondaries are strongly being irradiated. 
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1 INTRODUCTION 

Irradiation of the secondaries in close binary systems can 
significantly change their appearance and even their inter- 
nal structure. These are particularly dramatic in the case 
of low-mass X-ray binaries, in which the luminosity of the 
compact object can be as high as lO'^^ergs"^ (i.e. close to 
the Eddington luminosity of a neutron star), and where the 
luminosity intercepted by the secondary may be several or- 
ders of magnitude larger than the internal luminosity. De- 
pending on the structure and chemical composition of the 
stellar atmosphere, a fraction of the incident radiation will 
be absorbed, thereby heating material in the region of pen- 

* E-mail: podsi@astro.ox.ac.uk 



etration. Some of the luminosity deposited in the secondary 
will be re-radiated at lower wavelengths, increasing the out- 
ward flux and raising the surface temperature. The rest will 
move to cooler regions of the surface via circulation currents. 
This leads to a change in the temperature and entropy pro- 
files of the envelope. The former will affect the emitted flux, 
altering the phase-dependent lightcurve. It will also cause a 
variation in the temperature gradient, which can modify the 
spectral line profiles and, in turn, the radial velocity curves 
determined from spectral lines. The change in the entropy 
in the layer where the irradiation flux is thermalised may in- 
fluence the long-term equilibrium structure of the star and 
can, under certain conditions, lead to signiflcant expansion 
(Podsiadlowski 1991). 

Many of these effects have been addressed by a number 
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of authors in the past (see, for example, reviews by Pera- 
iah 1982; Kopal 1988; van Paradijs 1998; and also Davey 
& Smith 1992; Martin & Davey 1995). The purpose of the 
present study is to investigate an associated, but often ne- 
glected effect of irradiation: the pressure force due to the 
momentum transfer of the incident photons. This is partic- 
ularly important for X-ray binaries in which the radiation 
pressure force acting on the secondary may be a substantial 
fraction of the gravitational force. The consequences of this 
can be quite dramatic, sometimes involving a large defor- 
mation of the secondary's surface and a possible shift in the 
position of the effective inner Lagrangian point Li. In ex- 
treme cases, mass loss from the secondary can occur through 
one of the outer Langrangian points, L2 or L3, rather than 
the inner one. Irradiation can then drastically change the 
dynamics of mass transfer and the nature of the binary in- 
teraction. 

Standard procedures for the modelling of photometric 
lightcurves of close binary systems and most models for 
binary interactions are based on the Roche model for the 
equipotential surfaces of the binary components. This as- 
sumes that both components are centrally concentrated and 
that their gravitational effects can be approximated as those 
of point masses. The Roche model has been extended by 
various authors to overcome some of the intrinsic approx- 
imations of circular orbits and synchronous rotation (e.g. 
Plavec 1958; Kruszewski 1966). In addition, tidal and ro- 
tational distortions, the reflection effect, limb- and gravity- 
darkening can all be treated in a reasonably physical manner 
(e.g. Kopal 1959; Al Naimiy, 1978.). However, the classical 
Roche model allows only for gravitational and centrifugal 
forces, and no provision is made for systems in which exter- 
nal radiation pressure acting on either component is impor- 
tant. 

In this paper we present a systematic investigation of 
the effects of external irradiation pressure on the standard 
model of binary stellar surfaces. In section § 2 we justify the 
approach and in § 3 we show how the Roche-lobe model can 
be modified to approximately include the effects of external 
irradiation. In § 4 we describe the numerical procedure used 
to construct these approximate irradiated equipotential sur- 
faces. Binary systems with and without accretion discs are 
considered, and the validity of our solutions for given pa- 
rameters will be examined. In § 5 and § 6, we apply these 
solutions to a number of physical phenomena and systems, 
with and without accretion discs, respectively. In particu- 
lar, we discuss applications to binary pulsars, systems with 
extreme mass ratios and low- and high-mass X-ray binaries. 



2 THE RADIATION-PRESSURE FORCE 
TERM 

Radiation pressure is caused by the interaction between elec- 
tromagnetic radiation and stellar matter. Its strength de- 
pends on the momentum transfer per photon absorbed or 
scattered in the irradiated photosphere and is generally a 
very complicated function of local opacity and optical depth. 

In the presence of a radiation source, the equation of hy- 
drostatic equilibrium relating the pressure P, the density p, 
and the potential $ for any point r(r, 6, (j)) in the secondary. 



will acquire an additional term corresponding to the force 
due to this external radiation pressure, 

VP = -pV$+ ^/(r)ni, (1) 
c 

where ni is a unit vector in the direction of the exter- 
nal radiation, k is the local opacity, c the speed of light, 
/ the irradiation flux and r the local optical depth in the 
secondary envelope with respect to the external radiation. 
(Here we have neglected inertial terms due to the rotation of 
the system, i.e. the Coriolis force and circulation terms.) For 
a given source luminosity, /(r) at a point below the atmo- 
sphere will depend both on the local physical and chemical 
conditions in the secondary's envelope and on the geometry 
of the system, i.e. the distance from the radiation source 
and the angle of incidence of the flux vector. The compli- 
cated functional dependence of this term renders a precise 
analytical representation impossible. To obtain an explicit 
expression for this term, we generally make the simplifying 
assumption that the X-ray flux decreases exponentially (see 
e.g. Podsiadlowski 1991, 1992), 



47r|r — ri| 



(2) 



where \r — ri| represents the distance from the radiation 
source, and Lx denotes its X-ray luminosity. 

However, in general the irradiation term in equation 
(^) does not represent a conservative force, i.e. a force that 
can be written as the gradient of a scalar potential, since it 
depends on the local conditions in the region of flux absorp- 
tion. Therefore the inclusion of radiation pressure in the 
force equation precludes the existence of an equipotential 
surface. This also implies that there exists no surface for 
which the tangential force component is identical zero ev- 
erywhere. An accurate determination of the irradiated sec- 
ondary surface would therefore require a full treatment of 
the surface motion, including the meridional circulation of 
stellar material due to rotation (Von Zeipel 1924a) and the 
additional circulatory currents driven by the external ra- 
diation. A full solution is beyond the scope of this paper, 
and the correct circulation patterns even for non-irradiated 
single stars are still subject to debate (Osaki 1982; Char- 
bonneau 1992; Maeder & Zahn 1998). 

The purpose of this paper is to investigate the effects 
of the external irradiation in the absence of circulatory cur- 
rents, i.e. we consider the extreme case in which tangen- 
tial forces on the irradiated surface are small. In effect, we 
will seek a solution for which the irradiated surface locally 
approximates an equipotential surface. Compared to a so- 
lution that fully accounts for irradiation-driven circulation, 
one may expect that the true surface of an irradiated star 
will lie somewhere between this extreme case and the stan- 
dard Roche potential. 



3 A MODIFIED ROCHE POTENTIAL 

Even though the surface of an irradiated star cannot be rep- 
resented by an exact equipotential surface, one can still mod- 
ify the classical Roche potential to account for the radiation- 
pressure term in an approximate way. 

To to this, we flrst rewrite equation (|l|), expressing the 
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balance of forces per unit volume at the secondary surface, 
using the approximation of equation 

GMi 



+ 



\r — ri\^ 

Kp Lx 



ni 



c 47r|r — ri| 



(3) 



Here the effective potential <l>i combines the gravitational 
potential of the secondary, treated as a point mass, with 
the non-inertial potential terms due to the rotation of the 
system (i.e. due to the centrifugal force). 

Depending on the irradiation spectrum and the local 
physical and chemical conditions in the secondary, the ex- 
ternal flux will penetrate to a certain depth below the pho- 
tosphere. In our treatment of the irradiated surface, we will 
consider a layer of sufficient depth that most external radia- 
tion has been absorbed. However, this depth should be small 
compared to the secondary's radius, in order that |r — ri 
can be considered constant in this region and k and p can 
be well represented by average values R and p. Momentum 
transfer can then be regarded as an 'on the spot' event. 

We now integrate equation in the direction of the 
flux vector ni, from the top of the secondary's atmosphere 
to some fixed surface layer of vertical optical depth ro,n 3> 
1. We denote the path length of penetration as so with a 
corresponding optical depth to X-rays as ro = ro,nSec7, 
where 7 is the angle between the surface normal and the flux 
direction n\. With the assumption that only the exponential 
term changes signiflcantly over this distance, the right-hand 
side of equation (M) becomes 



GAfi 
r — ri I 



As 



Kp 



c An\r — ri| 



e ^ds. 



(4) 



The physical depth of penetration and the optical depth 
are related by dr = RpAs and therefore ro = Rpso- So for 
TO ^ 1, the right-hand side approximates to 

Lx 



GMi 
r — ri r 



^ ^P'\ 

r — r 

GMi 

= ~P'\ 1 

r - ri 



4:nc\r — ri P 
1 - 



AnG Micro 
ftK cos-y 



(5) 



where ft is the total integrated flux and g the gravitational 
acceleration due to the compact radiation source. Since both 
of these quantities have the same \r — ri\ dependence, the 
term in brackets in equation ^ is effectively constant for a 
given flux direction. Thus, at a given depth parallel to the 
top of the atmosphere, we may combine the forces due to 
gravity and radiation pressure from the compact object as 
a 'reduced' gravitational force 



Fgia,v — Fgra.v ~ -Frad — (1 — S)Fg 



(6) 



where S = constant x cos 7. The factor S, expressing the ratio 
of radiation to gravitational forces, therefore depends only 
on the cosine of the angle between the surface normal and 
the flux vector. However, it is this dependence that renders 
the reduced force non-central. 

We now extend the concept of a reduced force to de- 
flne a reduced potential. Both gravitational and radiative 
forces due the X-ray source can then be characterised as 



— {1 — 5) J F.dr = (1 — (5)0grav, and the standard Roche ex- 
pression for the binary potential may be modified directly 
by the inclusion of this (1 — S) factor. This formulation is 
strictly only valid in the case where gravitation and radia- 
tion pressure are acting in exactly opposite directions, and 
the resulting equipotential surface lies normal to both. The 
cos 7 dependence of 5 implies that this condition is not ful- 
filled, i.e. that a fictitious force results from the inclusion of 
a variable 5 factor in the total potential (see § 4.1). 

A number of authors have previously attempted to de- 
rive equipotential surfaces based on a reduced gravitational 
potential (e.g. Schuerman 1972; Kondo & McCluskey 1976; 
Vanbeveren 1977; Zhou & Leung 1988). However, these at- 
tempts did not fully account for the true physical and ge- 
ometric conditions of irradiated systems, resulting in unre- 
alistic solutions. We base our approach on the method de- 
veloped by Drechsel et al. (1995), who treated the geometry 
in a far more sophisticated way, although we use a different 
treatment for the inner Lagrangian point and extend their 
model to include the effects of an accretion disc, if present. 



4 DESCRIPTION OF THE MODEL 

Following Drechsel et al. (1995), we use a modified Roche 
potential that approximately accounts for the effects of ex- 
ternal irradiation to calculate the surface of the secondary, 
assumed to coincide with an equipotential surface. Since the 
external pressure force cannot be represented analytically, 
the equipotential surface must be calculated numerically. In 
the case of an accreting binary system, the model must also 
take into account the presence of an accretion disc, which 
may shadow a significant portion of the secondary and most 
importantly the inner Lagrangian (Li) point. In addition, 
when applying these modified surfaces to the calculations 
of photometric hghtcurves or radial velocity curves, the ef- 
fects of gravity- and limb-darkening, which cause significant 
perturbations to the surface temperature distribution, must 
also be included. This also requires a precise numerical de- 
termination of the gradients of the potential on the stellar 
surface. 

We consider the compact object as an isotropic point 
source of radiation, a fraction of which is incident on the sur- 
face of the secondary. The dimensionless ratio of radiation- 
pressure force to the gravitational force due to the compact 
companion at a given position on the secondary is defined 



5{r,e,4>) = 



(7) 



In the case where the irradiation flux is totally absorbed 
in the secondary's surface layer, which will be largely jus- 
tified in the case of irradiation by hard X-rays (see §4.4), 
these forces will be co-directional. We may then use the ap- 
proximations of a reduced gravitational force and a reduced 
potential introduced in the previous section. We take the 
centre of the secondary star as the origin of spherical polar 
coordinates, with the compact object along the polar axis 
of the coordinate system (i.e. in the direction 9 — 0) and 
the pole of the secondary in the direction 9 = 90°, = 0. 
The value of delta varies as the angle 7 between the surface 
normal and flux vector increases from its minimum value 
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obtained when the flux vector hes along the Une of centres, 
up to 90° at the horizon, causing a respective decrease in 
the radiation-pressure force proportional to the cos 7 term 
in equation (5). 

In the Roche model, the stellar surfaces are obtained as 
closed equipotential surfaces for given values of the poten- 
tial and mass ratios (Kopal 1959). When radiation pressure 
is included as a reduced gravitational potential, the modi- 
fied total potential for a surface element of the secondary 
star, assuming synchronous rotation, takes the dimension- 
less form 



1 1-S(r,0,(t>) 
= -+q , ' 
r Vl - 2Ar + 



q\r 



(8) 



where A = cos 9 and v = sin cos cj). 

The variation of 5 across the surface of the secondary 
has to be evaluated numerically. Once specified, the stellar 
surfaces can be determined point by point using a standard 
Newton- Raphson scheme (see e.g. Press et al. 1992). For 
given values of q and the critical potential f2o, the radius 
ro{9,(j)) for any direction {9,(j)) is obtained iteratively from 
the condition that the function g, defined as 



g{9,(j),r,no,q) ^Qo- ^{9, (j},r,q), 



(9) 



is identical zero. 

Writing 5 = 5max cos 7, we may relate 5niax to the stel- 
lar mass and luminosity of the conrpanion, and to the mean 
absorption coefficient R (see equ. H). This will generally be 
a complicated function, varying significantly across the sur- 
face of the secondary. In the present, exploratory investi- 
gation, we shall therefore make the simplifying assumption 
that (5max can be treated as an independent parameter. This 
has the added advantage of only introducing one new pa- 
rameter into the expression for the modified potential. This 
is desirable given the already large number of parameters 
necessary to specify the physical and geometrical structure 
of close binaries even for the classical Roche model. In §4.4, 
we shall make some rough estimates of the magnitude of 
(Jmax for typical parameters found in X-ray binaries 



4.1 Validity of the model 

Before we proceed to apply this model, we need to discuss its 
limitations. As already stressed in §2, the external radiation 
pressure force is not a conservative force and thus cannot 
be expressed precisely as the gradient of a potential. By 
assuming that the surface of the secondary coincides with 
one of these modified equipotential surfaces, we therefore 
make an error in the force equation (equ. 1) by effectively 
introducing a fictitious force. In addition, the assumption 
of hydrostatic equilibrium ignores the effects of irradiation- 
driven circulation. 



4-1.1 The fictitious force 

In our formulation, the external radiation pressure force is 
introduced as a reduction in the gravity term, (1 — 5)i^grav 
(as in equ. H). The modified potential due to the combined 
forces of radiation and gravity from the compact object is 
then taken to be —(1 — S) J F.dr, as in the second term 



of equation (g). However, this clearly introduces a fictitious 
and unphysical force due to the variation of 5 given by 



(10) 



In order to see the effect of this term on the determina- 
tion of the surface geometry, it is helpful to visualize that 
the construction of the equipotential surface is equivalent to 
piecing together surface elements that are perpendicular to 
the net local force. This provides a unique solution as long 
as there are a sufficient number of boundary conditions that 
allow a complete covering of the surface by these surface el- 
ements. In the present problem this is easily fulfilled since 
many surface points are already fixed, irrespective of irra- 
diation - for example, those in the shadowed region, once 
the potential at the inner Lagrangian point is specified. The 
forces derived from the modified potential then provide a 
good representation of the actual forces if, for all points, the 
fictitious force component does not significantly perturb the 
surface normal vector parallel to the actual force, — (Vfi)^, 
given by the gradient of the potential with 6 kept constant. 

To quantify this condition, we define the 'true' irradi- 
ated surface by the constraint 



{Vn)s A no = 0, 



(11) 



where no is the unit vector normal to the surface. By using 
a reduced potential, as given by equation (^, the surface 
actually calculated is given by 



[-{Vn)s + Fs] A [no + Sno] = 0. 



(12) 



The quantity Sno represents the deviation from the true sur- 
face normal and hence should be small (^ 1) for the model 
to represent the surface geometry faithfully. Expanding this 
expression and using equation ([ll[), we obtain 



[-(V!:!)i A Sho] + [Fs A ni] = 0, 



(13) 



where ni = no + Sho is the calculated surface normal in the 
model. But — (Vf2){ = |(Vf2)a|no, and hence we may write 



no A Snol — 



\Fs A ni 
\{VQ)s\ 



(14) 



It follows from the definition of a unit vector that 
I no A 5no| « \Srio\ (correct to second order in Sno). Thus, 
the requirement that |5no| <C 1 may be expressed as 

When this inequality is satisfied, the inclusion of the effects 
of external radiation pressure in the form of a modified po- 
tential only introduces a small error in the determination of 
the surface geometry. 



4.1.2 The role of circulation 

In our present model we do not consider the effects of 
irradiation-induced circulation in the outer layers of the sec- 
ondary (see e.g. Kippenhahn & Thomas 1979). Indeed, in 
general there is no physical reason why the surface of the 
secondary should be an equipotential surface, even approx- 
imately, if circulation terms are important. If a Roche-lobe 
filling secondary is suddenly being irradiated, its surface 
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would not immediately (i.e. on a dynamical timescale) ad- 
just to one of our modified equipotential surfaces, since the 
pressure inside the star, a few pressure scale heights below 
the surface, is always much larger than the external irradi- 
ation pressure. Instead, the main effect of a sudden turn-on 
of the external irradiation is to drive circulation that would 
move matter in the surface layer away from the irradiated 
side. However, this circulation can only be maintained as 
long as matter is resupplied from below. This will generally 
not be possible for an arbitrary surface geometry. Hence the 
surface will readjust until a closed steady-state circulation 
system is established (note, however, that even the existence 
of a steady state is by no means guaranteed). This demon- 
strates that, for a realistic treatment, the surface geometry 
of an irradiated secondary has to be determined simulta- 
neously with the circulation. The inclusion of circulation is 
beyond the scope of the present paper, but we have already 
initiated a detailed study of its consequences. By assuming 
that the surface is given by a modified equipotential, we ef- 
fectively minimize the effects of circulation. In this sense, 
these surfaces represent an extreme limit, and we may ex- 
pect that the true surface of an irradiated secondary lies 
somewhere between this limit and the other extreme, the 
unperturbed Roche potential. 



4.2 Evaluation of 5{r, 6, (f>) 

In our calculations we divide the surface of the secondary 
into an array of grid points spaced equally in 9 and (j). The 
value of S{r, 6, (p) is then determined for all grid points, where 
we need to take into account shadowing by an accretion disc, 
if present. In all cases, we assume that the compact object 
is a point source of radiation. 



4-2.1 Systems without accretion discs 

In systems without an accretion disc, there is no shielding of 
the secondary, and S has its maximum value along the line 
of centres, at the point 9 = (j> = 0. At all other points, S is 
calculated according to 



5{r,9,^)^ 



(5max cos7(r, 




if COS 7 > 
otherwise, 



(16) 



where j{r, 9, (f>) is the angle subtended by the surface nor- 
mal and the vector directed from the point (r, 9, cf>) to the 
source of the radiation. The constant (5max is an input pa- 
rameter, which defines the maximum value of 5 obtained in 
the absence of an accretion disc when the flux vector is in- 
cident normally on the irradiated surface. (Note, however, 
that there may not be any region of the surface on which 
the flux vector is incident normally; see §5.1.1.) 



4-2.2 Systems with Accretion Discs 

Low-mass X-ray binaries (LMXBs) in general and many 
high-mass X-ray binaries (HMXB) are believed to contain 
some form of an accretion disc. A disc will cast a shadow 
on a significant portion of the secondary's surface and will 
generally shield the L\ point from external irradiation. 

To model the effect of this shadowing, we assume that 



the accretion disc is axially symmetric and has an in- 
ner opaque region and an outer region which transforms 
smoothly from opaque to transparent at the outermost edge. 

We define the opening (half) angle of the disc as a and 
assume that the disc is completely opaque within an angle 
/3i of the equatorial plane (/3i < a), where these angles are 
defined with respect to the compact object. 

In Cartesian coordinates, where the compact object lies 
along the x-axis and where the axis of orbital rotation and 
the axis of symmetry of the accretion disc are parallel to the 
z-axis, P is given by 



sin/3 = \z\/d, 



(17) 



where d denotes the distance from a point on the secondary 
{x,y,z) to the compact object, given by 



d — ((a — x) -I- y + z 



2.1/2 



(18) 



and where a is the orbital separation. 

To model the partially transparent region of the disc, 
we define a composite function T{r,9,(j))-- 




if /? < /3i 

if /3i < /3 < Q 
ff /3 > a. 



(19) 



We also tested a linear transparency function in the partially 
transparent region, but found that a smooth sinusoidal de- 
pendence improves the rate of convergence in the numerical 
iteration process. 

Using this penetration factor, values of 5{r, 9, </>) are de- 
termined over the surface of the secondary in a similar man- 
ner to equation (lid). 



S{r,9,4>) = 



<5maxTcos7(r, 




if COS 7 > 
otherwise, 



(20) 



where the constant 5max is defined as before. 



4.3 The convergence procedure 

We now have a recipe for calculating S{r, 9, (j)) in terms 
of 7(r, 9, (j)) over the surface of the secondary. The angle 
7(r, (j)) is the angle subtended by the surface normal and 
the vector from the point (r, 9, (j)) to the radiation source. Its 
evaluation requires the knowledge of the 3-dimensional sur- 
face of the secondary under the influence of radiation pres- 
sure, which is a priori unknown and can only be determined 
iteratively. 

Following Drechsel et al. (1995), we model the surface 
of the secondary in the first iteration as a triaxial ellipsoid, 
defined by 



+ 



pole 



(21) 



with semi-axes rpdnt, ^gido, and rpoie identical with the 
Roche radii in the directions {9 — 0), {9 — 90°, — 90°), 
and {9 — 90°, (f} — 0), respectively. The radii rsidc and rpoic 
are taken to be equal to the classical Roche radii since these 
points lie beyond the radiation horizon. In the presence of an 
accretion disc, we assume that the inner Lagrangian point 
is shielded from radiation - hence rpoint will also be equal to 
the classical Roche radius. In systems which do not contain 
a disc, the position of the Li point will be modified. The 
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value of 5 at this point must simultaneously satisfy the re- 
quirement that the total force along the x-axis vanishes, i.e. 
dQ/dx = (where f2 is given by equation and must also 
satisfy equation (p^, which relates 5 to the surface nor- 
mal. Note that, as in the case of the unperturbed Roche 
lobe, the surface normal at the Li point is strictly unde- 
fined. We therefore consider a point infinitesimally close to 
the Li point and obtain an approximate surface normal by 
extrapolation. The value of rpoint is then fixed and forms 
one boundary condition for the iteration scheme. 

In each iteration, we need the vector from each surface 
grid point to the radiation source and the vector defining the 
surface normal. The former is readily calculated. However, 
a correct calculation of the surface normal is a crucial factor 
for a numerically stable iteration process, since this affects 
the variation of S with 9 over the irradiated region. For a cor- 
rectly converging solution, 5 must fall to zero continuously 
(though not necessarily smoothly) at the interface points, 
i.e. those points which separate the irradiated and shadowed 
regions. This forms a second boundary condition. We also 
require a continuous variation of S over the region shadowed 
by the accretion disc (this is ensured by our choice of the 
functional form of the transparency function in equ. 19) . 

Thus, for any point on the irradiated ellipsoid, the sur- 
face normal and flux vectors are determined. A value for 7 
can then be calculated and hence the local value of 5. Sub- 
stituting these S values into the modifled potential equation 
(^) then allows the determination of improved radius values 
for all points on the surface by application of the Newton- 
Raphson scheme. 



4.4 Physical estimates of (5max 

So far we have considered Jmax as a free parameter. However, 
it is straightforward to estimate the value of (5niax for actual 
binary systems. 

Consider an irradiated secondary. The radiation pres- 
sure acting on a unit surface element of the secondary is 
given by 



^max — 



-frad — 



/[/ COS 7 du dtu, 



(22) 



where c is the velocity of light, 7 is the angle between the 
surface normal and the incident radiation, du is the open- 
ing solid angle of the radiation cone, and the intensity 
in the frequency interval dv around v. Assuming a plane 
parallel atmosphere and negligible scattering of the incident 
radiation, we obtain for the radiation force per unit mass 



IdR, 



p dr 



Ki^/jycos 'ydh'duj, 



(23) 



where p is the mass density and k,^ the monochromatic ab- 
sorption coefficient for radiation absorbed by the secondary. 
By taking an integrated mean opacity, R, we may take this 
quantity out of the integral and obtain Jiad = ftii/c, where 
ft is the total integrated flux. The gravitational force due 
to the compact object per unit mass of the secondary is 
simply Fgrav = GM/d? = g, where M and d are the mass 
of the compact object and its distance from a point on the 
irradiated component. The maximum ratio of pressure to 
gravitational forces, assuming normally incident radiation 
(7 = 0), is therefore given by 



eg ' 



(24) 



which has the same functional dependence as in equation 
(P). We may also relate ^max to the Eddington luminosity 
of the secondary, I/Edd. Writing the flux as ft = L/iird'^, 
where L is the irradiating luminosity, 5max becomes 

Lk Lk 1 L 



And^cg AttcGM q I/Edd ' 



(25) 



where m is the mass of the secondary, and q — M/m is the 
mass ratio. 

In general, for a given frequency of radiation, the opac- 
ity will be a very sensitive function of the local tempera- 
ture, density and chemical composition, and can vary by or- 
ders of magnitude across the surface of the secondary. If the 
material on the surface of the secondary consisted of fully 
ionized hydrogen, electron scattering would be the domi- 
nant opacity source, i.e. R — aT/nip ~ 0.4 cm^g~^, where 
(XT ~ 6.7 X 10~^^ is the Thomson cross section. 

For systems where the luminosity in X-rays is close to 
the Eddington luminosity or where R is much larger than 
the Thomson opacity, the ratio of radiative to gravitational 
forces may easily be of order unity, implying a significant 
deformation of the classical Roche lobe. Indeed, it is possible 
that 5max exceeds unity. In this case, there would be no inner 
Lagrangian point. However, this is unlikely in practise as 
such systems are very likely to contain a thick accretion disc 
which will shield much of the inner face of the secondary. 



5 THE EFFECTS OF EXTERNAL RADIATION 
PRESSURE ON THE BINARY STRUCTURE 

External radiation pressure changes the structure and ap- 
pearance of the secondary compared to the case of Roche- 
lobe overflow without irradiation; the most important ef- 
fects are: (1) the distortion of the secondary's surface, and 
the associated change in the surface temperature distribu- 
tion which affects the modelling of lightcurves and radial 
velocity curves (see § 6); and (2) a possible change of the 
critical conflguration for systems which do not contain ac- 
cretion discs (or contain discs which are tilted with respect 
to the binary plane), where the Li point is directly irra- 
diated, which affects mass transfer and mass loss from the 
system. We will discuss these various possibilities in detail 
in this and the subsequent section. 



5.1 The binary configuration 

Let us first consider the effects of radiation pressure in sys- 
tems without accretion discs. This case has possible applica- 
tions to binary pulsar systems, stars orbiting supermassive 
black holes and, in some cases, high-mass X-ray binaries. 

To illustrate how external irradiation can change the 
binary configuration, we considered a binary with a fixed 
mass ratio q = M1/M2 — 2, but with different values for 
5max ~ 0.0, 0.15, 0.30 and 0.50 and calculated the critical 
modified potential, as described above, such that the sec- 
ondary fills its maximum possible volume. In Figure 1 we 
present meridional sections in the x-z plane of the modified 
critical potential surfaces for this sequence of models, where 
the same scale has been used in all panels. 
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(a) 6 = 0.0 




+ 



(b) 6 = 0.15 




(c) 6 = 0.30 




(d) 6=0.50 



-0.5 



0,5 1 -0.5 

X— coordinate (units of binary separation). 
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Figure 1. A sequence of modified 'Roche' lobes including the effects of external irradiation for a binary system with mass ratio 
q = M1/M2 = 2. The compact companion, of mass Mi, is represented as a single point at coordinates (1,0). The values of (5max, the 
maximum ratio of the external radiation force to the gravitational force due to the compact star at the secondary surface, are set to (a) 
0.0 (the standard Roche lobe), (b) 0.15, (c) 0.30, and (d) 0.50, respectively. 



Figure la shows the unperturbed Roche lobe without 
external irradiation pressure. In Figure lb the external irra- 
diation force has been increased to 15 per cent of the grav- 
itational force. As a consequence, the modified inner La- 
grangian point Li moves towards the compact object (since 
its effective gravity has been reduced), and the volume of 
the critical potential lobe increases. The critical potential, 
as defined by equation (^, decreases in magnitude as S is 
increased, until it reaches the same value as the potential 
of the outer Langrangian point, L2, which is unaffected by 
irradiation as it always lies in the shadowed region of the 
secondary. Once the outer Lagrangian point has a higher 
potential than the inner Lagrangian point, it takes over the 
role of defining the potential of the critical lobe, and the 
critical potential then passes through L2 rather than Li. 
This change from an inner to an outer critical configuration 
implies that mass loss now takes place through the outer 
rather than the inner critical point, most likely leading to 
the formation of a circumbinary disc and mass loss from the 
system rather than mass transfer to the compact object. 
For q — 2, this transition occurs at (5max = 0.202. In Figure 
Ic, where Jmax = 0.30, the switch of configuration has al- 
ready occurred. Increasing 5 further no longer changes the 
position of L2 (since its potential is unaffected by external 
radiation), but the shape of the critical potential is further 
deformed and its enclosed volume decreases (see Figure Id). 



5.1.1 The conical inner Lagrangian point 

As can be seen if Figure 1, for Jmax < (5crit, where Sait is 
the critical value at which the system configuration change 
occurs, the region around the irradiated Li point has a con- 
ical geometry (as for the standard Roche lobe). Hence the 
normal vector at the Li point is undefined, and its position 
must be determined from the asymptotic value obtained in 
its neighbourhood (i.e. in the limit of small 6, the angle 
measured from the line-of-centres, or x-axis) . Because of the 
centrifugal term in the total potential (equ. |^ there is no 
axial symmetry about the x-axis. The asymptotic value of 
the Li normal vector therefore depends on the plane con- 
taining the X-axis being considered, i.e. it depends on the 
angle (p. The critical equipotential surface is defined as the 
surface with the lowest potential for which any matter can 
leave the secondary's gravitational field. It is easy to show 
that this is the case for = 0; this implies that an expanding 
irradiated secondary star will first lose mass through the in- 
ner Lagrangian point in the x-z plane. (In practice however, 
if mass loss takes place, the overflowing material is likely 
to be optically thick and will therefore shield the Li point 
from radiation, thus obscuring the subtle geometrical details 
described above.) 

The conical nature of the irradiated Li point also im- 
plies that the maximum value of the ratio S, obtained by 
extrapolation, is less than (5max. At no point on the irradi- 
ated stellar surface is the fiux vector at normal incidence. 
This is in contrast to the model described by Drechsel et 
al. (1995) in which the Li point is apparently calculated 
assuming perpendicular incidence of the external radiation 
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flux, such that 5 = (5max at this point. Such a 'flat' L\ point 
also affects the critical value, iJcrit, at which the configura- 
tion change occurs. In the Drechsel model it is lower for a 
given mass ratio q. For example, for g = 2 it would occur at 
5crit = 0.167 for the assumption of a flat L\ point instead of 
5crit = 0.202, as obtained above. 

To show that a flat geometry near the L\ point is not ap- 
propriate, one can examine the factor -q defined in §4.1.1 (see 
equ. 15), which compares the fictitious force acting within 
the equipotential surface (resulting from the equipotential 
approximation) to the 'true' force directed approximately 
normally to this surface. We find that r] is at least an order 
of magnitude less than (5max at all points; this confirms that 
this modified equipotential surface provides a good repre- 
sentation of the 'true' forces. On the other hand, for a fiat 
Li geometry, r) is always of order unity in the neighbour- 
hood of the L\ point. This implies that the true and ficti- 
tious forces are of comparable magnitude and proves that 
the corresponding equipotential surface provides a poor ap- 
proximation for an irradiated surface. 

(Note that both models are consistent for outer critical 
configurations set by the L2 potential. In this case, the irra- 
diated surface at = is flat and 5 acquires its maximum 
value at this point, 5max.) 



5.2 Implications for binary pulsar systems 

As shown in § 5.1, if the irradiation force exceeds a critical 
value of the gravitational force, ^max, the critical potential 
configuration changes from an inner to an outer configura- 
tion. However, this can only apply to systems in which the 
inner Lagrangian point is directly irradiated and hence can 
never be the case for systems like low-mass X-ray binaries 
where the accretion disc effectively screens the L\ point from 
the X-ray irradiation. In such systems, the main consequence 
of irradiation will be a deformation of the surface (see § 6). 

On the other hand, binary systems containing millisec- 
ond pulsars probably do not contain accretion discs. Hence 
pulsar radiation (in the form of electromagnetic waves, par- 
ticle flux, etc.) can play the role of X-rays and can poten- 
tially cause a conflguration change. Mass will then be lost 
from the outer Lagrangian point, which may then leave the 
system entirely. Continued irradiation may eventually lead 
to the total evaporation of the secondary. This is the main 
mechanism, originally proposed by Ruderman, Shaham & 
Tavani (1989), by which pulsars may be able to destroy their 
companions leaving a single millisecond pulsar. The "black- 
widow" pulsar, PSR 19574-20 (Fruchter et al. 1988), and the 
binary pulsar PSR 1718-1-19 (Lyne et al. 1990; McCormick 
& Frank 1993) are believed to be examples where the pulsar 
radiation is in the process of evaporating the companion. 

In the case of PSR 1957-1-20, the mass function and the 
radio eclipse suggest a secondary mass near 0.022 M© (for 
an assumed pulsar mass of IAMq; e.g. Fruchter et al. 1988). 
This implies a mass ratio of around 60 and a corresponding 
value for Sc-cn of approximately 0.01 (see Table 1). One can 
also estimate the pulsar luminosity from the observed spin- 
down rate to be ~ 401/0 (e.g. Fruchter et al. 1990). For the 
irradiated component we have I/Edd — (10^^ k~*) ergs~^, 
and therefore 5max = L/gl/Edd = 0.0025 k. So 5max will ex- 
ceed Jcrit if the mean photospheric opacity on the irradiated 
side is greater than ~ 4 cm^ g~*, a value that may not be 



Table 1. Values of <5criti 
qucncc of mass ratios 



the product 



and f2c 



for a se- 



Mass Ratio, q 






^crit 


0.5 


0.459 


0.230 


2.408 


1.0 


0.318 


0.318 


3.207 


1.5 


0.246 


0.369 


3.949 


2.0 


0.202 


0.405 


4.655 


2.5 


0.172 


0.431 


5.335 


10 


0.0564 


0.564 


14.44 


10^ 


6.89E-3 


0.589 


109.7 


10=* 


7.48E-4 


0.748 


1.02E-f3 


10'* 


7.76E-5 


0.776 


l.OOE+4 


10^ 


7.89E-6 


0.789 


l.OOE+5 


10** 


7.94E-7 


0.794 


l.OOE+6 



unreasonable (it depends on the uncertain thermodynamic 
parameters in the irradiated photosphere of the secondary). 
In this case, we may expect a configuration change where 
mass loss occurs from the outer Lagrangian point, leading 
to the formation of an excretion disc rather than an accre- 
tion disc. 



5.3 Application to binary systems with extreme 
mass ratios 

Table 1 shows the values of i5crit, the product g^crit, and the 
critical potential ficrit of the secondary, as defined by equa- 
tion (^), for a sequence of mass ratios. For higher mass ratios, 
the critical value of (Jmax decreases, and the magnitude of 
the critical potential increases, corresponding to a smaller 
Roche volume. Note that the product g^crit converges to- 
wards a constant value for increasing q. In conjunction with 
equation (^), this implies that in the asymptotic limit of 
high mass ratio, the critical luminosity scales directly with 
the Eddington luminosity of the irradiated object. 



Lc-cit — g^crit ^Edd 

oc LEdd, for g > 1. 



(26) 



A related formula was derived by Podsiadlowski & Rees 
(1994), where they expressed the difference in potential be- 
tween the inner and outer Lagrangian points in terms of the 
ratio of irradiating to Eddington luminosities. Their equa- 
tion (7) (without a spurious factor of (m/3M)*/^) reads 



Gm 



(27) 



where 7 — L/L^dd, a is the binary separation, and m the 
mass of the irradiated star. However, this equation contains 
the assumption that the irradiating flux is perpendicular to 
the irradiated surface at Li, /5li ~ 1. As we have 

shown above, this is not a good approximation. Based on 
our calculations, performed for a range of values of q and 
5max, we find this ratio to be approximately 1.20. Hence the 
Lagrangian points have the same potential, i.e. A?7 = 0, for 
7 — 1.20 X |. Equation (|2q) therefore becomes 



0.80I/Edd, 



(28) 
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in agreement with table 1. 

One interesting application of this case are low-mass 
stars orbiting supermassive black holes in active galactic nu- 
clei (AGN) which are being irradiated by the central black 
hole. Such systems can be considered binaries with extreme 
mass ratios, typically of order 10^ to 10*, and some highly 
unusual properties (see e.g. Podsiadlowski & Rees 1994). 

There is a general current consensus that AGN are pow- 
ered by the accretion of matter onto the central black hole 
(see e.g. the reviews by Rees 1984; Frank, King & Raine 
1992). Perhaps the strongest evidence for this picture de- 
rives from studies of X-ray emission, in which the combina- 
tion of large luminosity (typically ~ 10% of the bolometric 
luminosity) and rapid variability imply an efficient energy 
generation mechanism confined within a small region (e.g. 
Fabian 1992). Although accretion is the most plausible ex- 
planation, it is not clear what the source of the fuel is. Tidal 
disruption of stars is one possibility or even Roche-lobe over- 
flow of a star that spirals towards the black hole due to 
angular momentum loss caused by gravitational radiation 
(see Podsiadlowski & Rees 1994 for discussion and further 
references). 

However, the X-ray luminosity in these systems will 
generally be far in excess of the Eddington luminosity of 
the stellar companion, ~ 10'^* ergs~^ for a star of order 
1 Mq . (In general, for luminous supermassive black holes, 
we may expect L/L^dd ^ 10^-10^). So it is clear from equa- 
tion (^8|) that only a tiny fraction of the central source's 
X-ray luminosity has to be emitted in the direction of the 
companion and reach the inner Lagrangian point for L to 
exceed Lcrit . Hence any Roche-lobe overflow from the stellar 
companion must occur via the outer Lagrangian {L2) point, 
leading to the formation of an excretion disc rather than 
an accretion disc. The tidal coupling between this excretion 
disc and the secondary may then provide a tidal barrier and 
actually prevent accretion of disc material onto the black 
hole (e.g. Artymowicz et al. 1991). 



5.4 Self-sustaining mass loss through the outer 
Lagrangian point 

In the previous sections, we showed how external irradia- 
tion pressure can change the configuration of the critical 
tidal lobe for the mass-losing component in a binary sys- 
tem from an inner to an outer critical configuration where 
mass loss occurs through the outer rather than the inner 
Lagrangian point. Since this matter will carry off an associ- 
ated amount of angular momentum that is larger than the 
systemic specific angular momentum, this leads to a net de- 
crease in the specific angular momentum of the system. Such 
mass loss can dramatically increase the mass-loss rate from 
the mass loser and sometimes even destabilize such systems 
completely. To estimate this more quantitatively, we give a 
simple analysis following along standard lines (see e.g. Rap- 
paport, Verbunt & Joss 1983). 

We consider a binary system whose components have 
masses Mi and M2, and circular Keplerian orbits of radius 
oi and 02, with a total separation a. The system has an 
orbital period P, and the total orbital angular momentum 
of such a configuration is given by 



J : 



Ga 



Ml + M2 



M1M2. 



(29) 



Now suppose that the secondary loses mass by an amount 
5M2, a fraction /3 of which is accreted by the primary while 
the remaining fraction leaves the system. In units of the 
orbital angular momentum of the secondary, the angular 
momentum lost from the system is therefore 

5J^a{l-[3)5M2al'^, (30) 

where the factor a depends on the details of the mass- 
loss mechanism. Using Kepler's Law and the relation a2 — 
Alia/ (Ml + M2), we can combine equations ( p9| ) and ( |30[ ) 
to obtain 

Ml 



J M2 M1 + M2 



(31) 



Taking the logarithm of equation (g9|) and differentiating 
it with respect to M2, and then equating the result with 
equation (pll), we get 



5 In a 



= (1-/3) 



2aMi -I- M2 



(32) 



(5 In M2 Ml -I- M2 

Since the orbital separation changes because of the mass 
loss, the critical tidal lobe will adjust accordingly. Writing 
the effective radius of the critical tidal lobe as 



Rl = af{q), 



(33) 



where q = M-1/M2 is the mass ratio and the function / de- 
pends on the exact shape of the critical lobe (for the stan- 
dard Roche lobe, see the approximation by Eggleton [1983]). 
Even in the case of the irradiation-modified critical lobe we 
can factorize Rl in this form as long as the geometry of 
the modified potential remains self-similar (for a constant 
irradiation luminosity) . 

As mass is lost through the outer Lagrangian point, 
the separation will generally shrink and so will the critical 
lobe. If the critical radius decreases more rapidly than the 
thermal equilibrium radius of the secondary, the latter will 
no longer be able to remain in thermal equilibrium, and mass 
transfer will then be self-sustaining, i.e. be driven by the 
thermal expansion of the secondary or, in the most extreme 
case, occur on a dynamical timescale (see e.g. Ritter 1996 
for a detailed review). This condition is best expressed in 
terms of the mass-radius exponent of the critical lobe and 
the equilibrium radius as 



where 

dlni? 



/dlni?\ , ^ /dlnii\ 

^^=idh^jL' ""^^-=ldh^j. 



(34) 



(35) 



define the mass-radius exponents for the critical tidal lobe 
and stars in thermal equilibrium, respectively. Using equa- 
tions d33) to (pa), we then obtain 



Cl = -2 + (1-/?) 



2aq 



1 



dlnf 
din g 



(36) 



Assuming now that all the mass is lost through the outer La- 
grangian point and none is accreted by the compact object, 
and that this matter carries off exactly the specific angu- 
lar momentum at the outer Lagrangian point, located at 
x-coordinate xl^i we can set (3 = and obtain a from 



10 S. Phillips and Ph. Podsiadlowski 



1 + 



\XL2 



where 02 = Mia/ {Mi + A'h)- Equation 
to become 



(37) 



then simplifies 



-2 + 



2aq + l din/ 
din q 



1 



(38) 



The function din //ding is generally a slowly varying 
function which will depend on the irradiation luminosity. In 
the case of the standard Roche formula (Eggleton 1983) , it 
is easy to show that is > 1 for q < 100. Since is gen- 
erally less than 1, this suggests that mass loss will generally 
be self-sustaining for q< 100, occurring either on a thermal 
or a dynamical timescale, if mass loss occurs through the 
outer Lagrangian point. For comparison, in the case of con- 
servative mass transfer through the inner Lagrangian point, 
the condition — 1 corresponds to a mass ratio q ~ 0.8. 

The change of the orbital separation is related to a 
change in orbital period by Kepler's Third Law which can 
be written in its differentiated form as 



3 a 
2a 



(1 - /3)M2 



(39) 



we obtain 



2(Mi + Ma) • 
Combining this with equation 

Porb ^ „M2 
Porb M2 ' 

where 

r_P f 2/3 + /?/3 + g[l- 
<1 \ 1 + </ 

(This treatment does not include any tidal interaction terms 
due to a difference in rotational and orbital angular frequen- 
cies.) In the case of non-conservative mass loss through the 
L2 point (/3 = 0), this reduces to 



a{l~l3)] 



(40) 



(41) 



-Porb 
-Porb 



2-^3g(l-a) 



M2 
M2' 



(42) 



relating the period change to the mass-loss rate from the 
secondary. 

In recent years, a number of systems have been observed 
to exhibit orbital period decays, for example, the massive X- 
ray binary SMC X-1 with Porb/^orb = (-3.36±0.02) x lO"'^ 
yr~^ (Levine et al. 1993) and the system Centaurus X-3 with 
Porb/Porb = (-1.8 ±0.1) X 10"^ yr"^ (Kelley et al. 1983). 
These period changes have traditionally been explained by 
a tidal instability, the Darwin instability, which causes the 
orbit to decay in the process of spinning up the secondary 
(Darwin 1879; Pringle 1974). 

However, as we shall demonstrate below, using Cen X- 
3 as an example, a relatively moderate outflow through L2 
could also account for the observed change and may help ex- 
plain some other observations that are not well understood 
at the present time. 



1983), this is about 50% greater than the flux required to 
produce a configuration change provided that the inner La- 
grangian point is directly irradiated. Assuming that mass 
loss occurs through the outer Lagrangian point, we can cal- 
culate a from equation (36) to be ~ 460, where XL2 was 
calculated numerically. Equation (^2|) then becomes 



(43) 



Porb „„M2 

— — ~ 60-T-. 

Porb M2 

This should be compared to the case where mass transfer is 
conservative through L2, where the factor in front of M2/M2 
would be ~ 20 instead of 60. Thus, to reproduce the ob- 
served rate of orbital decay (Porb/Porb ~ —1.8 x 10~^ yr~^), 
a mass loss rate of around 6.0 x 10~^ M0 yr~^ through the 
outer Lagrangian point would be required. This is consis- 
tent with the initial rates expected for overflow through the 
Li point in HMXBs during the early phase of mass transfer 
(e.g. Lamers, van den Heuvel & Petterson 1976). In this case, 
however, one expects the formation of an excretion disc em- 
anating from the L2 point (see e.g. Shu, Lubow & Anderson 
1979), from which emission or absorption lines may be ob- 
served depending on the orientation of the observer. Indeed 
strong evidence for material extending beyond the limiting 
Roche lobe surface of the secondary has been reported by 
Mauder (1975) and Clark et al. (1988), although the true ori- 
gin of this material remains unclear. If indeed the secondary 
loses mass through the outer Langrangian point, an efficient 
mechanism would be required to fuel the accretion-powered 
X-ray source, most likely in the form of a stellar wind, quite 
likely enhanced by X-ray irradiation (Day & Stevens 1993). 

Figure 2 shows approximate equipotential surfaces for 
Centaurus X-3 for a mass ratio q = 0.05. The dotted curve 
represents the unperturbed Roche lobe, with the compact 
X-ray source at coordinates (1,0). The dotted inner region 
represents the deformed secondary when X-ray irradiation 
(assumed to be Eddington limited) is included, and the Li 
point is shadowed by a thin accretion disc (not shown) . The 
cross-section of the irradiated Roche lobe is significantly re- 
duced by the incident radiation, causing a bottleneck near 
the Li point and strongly reducing the effective cross section 
for a stream emanating from Li (see Lubow & Shu 1975). 

If the accretion disc is either tilted or absent, allowing 
direct irradiation of the Li point, the secondary can expand 
to fill the outer equipotential surface (solid line) . Matter will 
then preferentially escape from the outer Lagrangian point 
forming an excretion disc (indicated schematically as shaded 
regions) . 

If Cen X-3 is losing mass through the outer Lagrangian 
point, this has several important consequences both for the 
evolutionary stage of the system and the modelling of the 
system. Although we use Cen X-3 as an example, the model 
could apply equally well apply to similar high-luminosity 
HMXBs, e.g. SMC X-1. We require only that the X-ray fiux 
at the Li point is greater than the critical flux; for systems 
radiating close to the Eddington limit, this will generally be 
the case. 



5.5 Application to Centaurus X-3 

Centaurus X-3 is a massive X-ray binary with an orbital 
period of 2.87 d, observed to radiate near the Eddington lu- 
minosity for a neutron star, ~ 2 x 10^® (White, Swank & 
Holt 1983). For a mass ratio of 0.05 (Rappaport & Joss 



5. 5. 1 The lifetime of the X-ray phase 

One consequence of mass loss through the outer Lagrangian 
point is that it may increase the lifetime of the X-ray-active 
phase of the system. High-luminosity HMXBs are commonly 
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Figure 2. Approximate equipotcntial surfaces for the high-mass X-ray binary system Centaurus X-3, assuming a mass ratio of 0.05. The 
dotted curve represents the unperturbed Roche lobe, with the compact X-ray source at coordinates (1,0). The dotted inner region shows 
the deformed geometry of the critical lobe when the effects of X-ray irradiation are included (for an Eddington-limited X-ray source), 
but the Li point is shadowed by a thin accretion disc (not shown). The thick solid curve shows the critical potential surface when the 
Li is irradiated directly with a luminosity of 2 X 10^® ergs~^. Matter will then preferentially escape from the outer Lagrangian point, 
forming an excretion disc, as indicated schematically by the shaded regions. 



assumed to be powered by Roche lobe overflow through the 
Li point (Lamers et al. 1976). However, this will rapidly 
lead to runaway mass transfer on a thermal time scale, 
which would result in accretion rates of order 10"* to 10~^ 
M0yr~^ (van den Heuvel & De Loore 1973). Such large 
rates would quickly and completely extinguish the X-ray 
source by absorption, suggesting a maximum X-ray lifetime 
of < 3x 10^ yr (Terman, Taam & Savage 1998). In the case of 
outflow through L2, all of this mass loss from the secondary 
will feed a circumbinary disc. The X-ray active lifetime will 
be determined by the rate at which the wind accretion rate 
on the neutron star increases. It could be comparable to the 
lifetime estimated for low-luminosity HMXBs, ~ lO^'-lO^ yr, 
which are thought to be predominantly wind-fed (Meurs & 
van den Heuvel 1989). 



5.5.2 The size of the X-ray eclipse 

Cen X-3 is an eclipsing system where the total X-ray eclipse 
lasts 0.488 d (Schreier et al. 1972), i.e. a fraction of 0.17 of 
the total period. Deflning the eclipse half-angle 9c as the 
orbital phase angle at which the occultation of the X-ray 
source is just complete, the data implies 6^ « 42°. This 
angle may be directly related to the radius of the secondary 
R2 according to 

^2 /, . 2 . 2 o \l/2 

— — (1 — sm zcos Oc) , (44) 

(Rappaport & Joss 1983) where i is the system inclination 
and a the component separation. If one assumes that the sec- 
ondary flUs its standard Roche lobe, one can then use Eggle- 
ton's (1983) expression for the Roche lobe radius to uniquely 
determine the mass ratio q. Since the mass function is ap- 
proximately 15 Mq (Schreier et al. 1972), this would yield a 



mass of the compact object of less than 0.7 M© (Davidson 
& Ostriker 1973), which is not consistent with current esti- 
mates of ~ 1.4 Mq (e.g. Thorsett et al. 1993) for a typical 
neutron star. Clark et al. (1988) attempted to fit the Cen 
X-3 data using a spherically symmetric coronal wind model 
which took into account the gradual eclipse ingresses and 
egresses. Their model predicted a more reasonable neutron 
star mass of 1.23 ± 0.60 Mq, but required a much lower 
(total) eclipse half-angle of 32.9°. 

A much simpler solution of this discrepancy is to in- 
crease the size of the eclipsing region, as occurs naturally 
if the secondary is filling it outer critical lobe (see Fig. 2). 
To demonstrate this, we consider an inclination angle in the 
probable range 70-80° (Nagase 1989) and assume that nei- 
ther the star itself nor any material associated with it can 
exceed the dimensions of the irradiated Roche lobe. Follow- 
ing the analysis of Schreier et al. (1972), we adopt 6c ~ 45° 
(although Wilson [1972] has suggested values as high as 53°). 
Equation ^4) then implies that the irradiated Roche lobe 
radius is in the range 0.72-0.75, in units of the orbital sep- 
aration. Using our model, rather than Eggleton's formula, 
we relate R2 to q by calculating the total volume of the ir- 
radiated lobe for a given mass ratio and equating this with 
|7ri?2. We obtain values for q in the range of 0.101-0.062, 
which corresponds to a far more plausible range for the neu- 
tron star mass of 1.26-1.92 Mq. One can expect that a more 
thorough investigation will yield tighter constraints on the 
possible masses, the inclination and the eclipse angle. 



5.5.3 The optical lightcurve 

A second observational problem associated with Cen X-3 is 
related to the modelling of the optical lightcurve of the sys- 
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tem (e.g. Mauder 1975; Tjemkes, Zuiderwijk & van Paradijs 
1986). Using a standard model for the ellipsoidal variations, 
Tjemkes et al. (1986) found that the depths of the minima 
at phases 0.0 and 0.5 were much greater than the typical 
model parameters predicted. They discounted the possibil- 
ity of X-ray heating since this would further decrease the 
depth at phase 0.5. However, Hutchings et al. (1979) pre- 
sented conflicting evidence. They found the spectral type 
of the companion star to vary with orbital phase, ranging 
from 09 near phase 0.0 to 06 near phase 0.5, which strongly 
suggested appreciable X-ray heating. Our model would ap- 
pear to broadly reconcile these results if we allow irradiation 
of the Li point, since the initial increase in the companion 
Roche lobe actually causes an increase in this depth due the 
change in geometry. Clearly, this effect should be included 
in accurate lightcurve modelling of the system. One should 
also note that strong irradiation need not necessarily result 
in a strong 'reflection' effect, since only a fraction of the inci- 
dent X-ray flux has to be thermalised, while circulation cur- 
rents will further distribute the heated material around the 
secondary (Kippenhahn & Thomas 1979; Schandl, Meyer- 
Hofmeister & Meyer 1997). 

5.5.4 Caveats 

In our analysis we made two important assumptions: (1) 
that the secondary is in (almost) synchronous rotation, such 
that the Roche lobe approximation may be applicable; and 
(2) that the Li point is directly irradiated, with an associ- 
ated change in the value of the critical potential. 

It is generally believed that any initial non-circularity 
in the binary orbit or frequency difference in the rates of or- 
bital and axial rotation would have been rapidly removed by 
tides exerted by the compact star on the secondary. The syn- 
chronisation timescale in X-ray binaries is generally of the 
same order of (or shorter than) the circularization time scale 
(Lecar, Wheeler & McKee 1976; Zahn 1975; Hut 1981). As 
Fabbiano & Schreier (1977) have shown that the secondary 
in Cen X-3, also known as Krzeminski's star, is in a highly 
circular orbit (e = 0.0008± 0.0001), it is safe to assume that 
secondary has synchronized its spin with the orbit. 

The second of these caveats is somewhat more problem- 
atic. In the simplest model of an X-ray binary, one assumes 
that matter is accreted onto the compact object through an 
accretion disc which forms in the orbital plane. The Li point 
is then shadowed from a point-like X-ray source positioned 
at the compact object. However, there are a several possible 
alternatives. 

First, even if an accretion disc is present, the Li point of 
the secondary may be directly irradiated by X-ray radiation 
coming from an extended Compton-heated thick accretion 
disc corona for which evidence has been found in the partial 
X-ray eclipses of several non-pulsating accretion-powered bi- 
nary X-ray sources by White & Holt (1982). 

Second, the disc plane may be tilted with respect to the 
binary plane and allow direct irradiation of the secondary 
(depending on the precession phase). Such tilted discs have 
variously been proposed to explain precessing discs in sys- 
tems like Her X-1, LMC X-4 and SS 433 (see e.g. Roberts 
1974; Petterson 1977; Priedhorsky & Holt 1987) and even 
Cen X-3 (Iping & Petterson 1990). As Pringle (1996) has 
shown, even initially planar discs are unstable to warping 



due to the external radiation pressure, implying that warp- 
ing is a generic phenomenon in centrally illuminated accre- 
tion discs (also see Maloney & Begelman 1997). 

Third, it is by no means clear that an accretion disc 
forms if mass accretion onto the compact object occurs 
through a stellar wind, since very little of the orbital angu- 
lar momentum is actually accreted with the wind material 
(see e.g. the studies by Ruffert 1994; 1997). No disc is ex- 
pected to form if the circularization radius, which depends 
on the specific angular momentum of the accreted material, 
is less than the Alfven radius, which defines the radius below 
which the neutron star's magnetic field will determine the 
dynamics of the accretion flow. 



6 BINARY SYSTEMS CONTAINING AN 
ACCRETION DISC 

In general, LMXBs and some HMXBs contain substantial 
accretion discs. These shield the inner Lagrangian point 
from X-ray irradiation, thereby maintaining the relative po- 
sitions and values of the Lagrangian points. However, radi- 
ation pressure may still cause significant deformation of the 
irradiated stellar surfaces, particularly in systems with high 
accretion rates. This can change the appearance of the op- 
tical lightcurves and radial velocity curves, thereby altering 
the derived system parameters. 

6.1 Observational effects of surface deformations 

To illustrate the effects of radiation pressure, we consider 
a binary system with physical parameters that may be ap- 
propriate for the LMXB Scorpius X-1 with an orbital pe- 
riod of 18.9 hr (Kallman, Boroson & Vrtilek 1998), where 
we take the masses of the neutron star and the secondary to 
be 1.4 M0 and IM0, respectively. The system then has an 
orbital separation of 4.8R0. The irradiating X-ray luminos- 
ity is assumed to be 2 x 10^* ergs~^, which is approximately 
the Eddington luminosity of the neutron star. The optical 
companion in this system is thought to be a Roche-lobe 
filling star near the end of or just beyond its main-sequence 
phase (Cowley & Crampton 1975). We take an effective (po- 
lar) temperature of IGOOOKQ, a gravity darkening coefficient 
of 0.08, and a constant limb-darkening coefficient of 0.35, 
appropriate to V band observations (see appendix). For the 
accretion disc, we assume a conservative opening (half) an- 
gle of 10° with a fully opaque region of 7°. These parameters 
imply a mass ratio of g = 1.4, and therefore a ratio of radi- 
ation to pressure forces Jmax ~ L/qL-^dd ~ 2.5k. Using the 
approximation described in §4.4, we will assume k = 0.4, 
implying (5max = 0.99. We emphasize that this model should 
not be considered very realistic, since it makes a number of 
simplifications, e.g. it assumes constant opacity across the 
irradiated surface and that all the energy is thermalized be- 
low the photosphere of the secondary (i.e. takes a X-ray 

^ This temperature is significantly higher than the eff'ective tem- 
perature of a 1 Mq star, but could be appropriate if the secondary 
were in a similar evolutionary phase as the secondary in Cyg X-2 
(Casares, Charles &; Kuulkers 1998) or if surface circulation were 
very efiicient in re-distributing the irradiation flux around the 
secondary 
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Figure 3. Deformation of the Roche Lobe for 5max = 0.99, 
q = 1.4. The outUiie of the standard Roche lobe is shown as 
a dashed curve, and tlie X-ray source is represented as a point 
with coordinates (1,0). The cross-section of the accretion disc is 
shown containing a heavily shaded opaque region and a more 
lightly shaded transition region. 



standard Roche Model 
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Figure 4. Comparison of V band optical liglitcurves due to stan- 
dard and modified Roche lobes, for the parameters given in the 
text. Relative magnitudes are calculated by subtracting the theo- 
retical magnitude at each phase, m, from the magnitude at phase 
0, mo. 



albedo of 1). The main purpose of the model is to demon- 
strate the magnitude of the effects of irradiation pressure on 
the derived system parameters. 

The results of the model using these parameters are 
shown in Figure 3, together with the outline of the Roche 
lobe due to the standard model. The compact object is 
again represented eis a single point with coordinates (1,0). 
The cross-section of the accretion disc is shown, containing 
a heavily shaded opaque region and a more lightly shaded 
transition region. Its effect as an X-ray shield may be readily 
inferred from the lack of any distortion around the Li point. 

The deformation of the Roche surface is very clear. The 
effect is greatest in the region just beyond the accretion 
disc shadow, and decreases as the Roche lobe horizon is ap- 
proached. It should also be noted that the surface equipoten- 
tials 'cuts' the interior equipotentials to a depth determined 
by the extent of the deformation (about 1/5 of the radius 
for the illustrated case). This may have implications for the 
circulation of energy along equipotentials. 

6.1.1 Effects of radiation pressure on optical lightcurves 

To demonstrate the effect of X-ray radiation pressure on 
optical lightcurves, we have calculated these (using the 
lightcurve code described in the appendix) for both a stan- 
dard Roche lobe model and a model with a distorted Roche 
lobe (Figure 3) for the parameters given above. In order to 
avoid complications due to eclipses of the secondary or the 
accretion disc, we chose an inclination of 40°. Figure 4 shows 
the resulting optical V band lightcurves calculated by sum- 
ming the flux over the surface of the modified Roche lobe 
and the standard Roche lobe, respectively. In both cases, 
the accretion disc is assumed to contribute a fixed amount 
of flux. For each curve, the relative magnitude is plotted 
against phase, where the former is calculated by subtract- 
ing the theoretical magnitude at each phase, m, from the 
magnitude at phase 0, mo. This vertical shift ensures that 
both curves pass through the origin, thereby highlighting 
their differences at phase 0.5. 

The greatest difference between the two models occurs 



at phase 0.5, for which the irradiated side of the secondary 
is most visible and the observed flux is a maximum. The flux 
due to the pressure-modified lobe is lower, since the distor- 
tion results in a decrease in the heated area perpendicular 
to the line-of- sight. Around phase 0, the curves differ much 
less, since the backside of the secondary is visible, and its 
geometry is not affected by radiation pressure. 

The amplitude of a lightcurve is defined as the differ- 
ence in magnitudes between phases and 0.5 and is gener- 
ally a sensitive function of the inclination of the system. It 
provides an empirical correlation between the observations 
and the system parameters (e.g. Bochkarev, Karitskaya & 
Shakura 1979). If the mass ratio q is known for a system, 
for example from radial velocity curves determined spectro- 
scopically, it is common to use the lightcurve amplitude to 
deduce the orbital inclination and with it the component 
masses, which arc proportional to sin i to the third power. 
Thus, a fractional error in the value of sini propagates to 
three times this error in the component masses: 

Am 3A(sini) 

^ ^— 45 

m s\m 

The top panel of Figure 5 shows the theoretical ampli- 
tudes expected for this system as a function of inclination, 
calculated for both the standard Roche model and the irradi- 
ated model. At an inclination of zero, the observer views the 
system face-on an no variation in the lightcurve is observed. 
The amplitude has a maximum at an inclination of 90° when 
the line-of-siglit is parallel to the orbital plane (of course, at 
large inclination eclipse effects will also strongly affect the 
lightcurve). As described earlier, the theoretical amplitude 
is systematically smaller for the irradiated model than in 
the standard model. Thus, for a given observed amplitude, 
different models will correspond to different inclinations and 
hence implied masses (through equ. 45). In the lower panel of 
Figure 5, we have plotted the quantity 3 A sini/ sini against 
inclination, where we assume that the true value of sin i 
refers to the irradiation model. This correction is significant 
(> 7 per cent) at all inclinations and rises beyond 30 per cent 
at the extremes. At an inclination of 40° to the line-of-sight, 



14 S. Phillips and Ph. Podsiadlowski 



\ 0.2 



0.1 



Standard Roche Model 
Irradiation Model 




inclination (degrees) 

Figure 5. Top panel: Theoretical V band amplitudes obtained 

for tlic standard Roche model and irradiation model as a function 
of inclination. Bottom panel: The fractional error in (sini)^ as 
a function of inclination. For a given amplitude, the irradiation 
model is assumed to provide the true inclination. 



the value of 3Asini/sini is approximately 0.26. implying 
a, ^ 25 per cent error in the component masses. At lower 
inclinations, the error is greater still. 

Systems where external radiation pressure can be ex- 
pected to be important are systems such as Scorpius X-1, 
Cygnus X-2 or the black hole candidate GRO J1655-40, in 
which the secondary is sufhciontly evolved so as to allow op- 
tical observations during outburst (Orosz & Bailyn 1997). 
The effect will also be of importance in HMXBs, whose op- 
tical output is dominated by the high-mass secondary star 
(for example, the super-Eddington X-ray binary, LMC X-4; 
see Heemskerk & van Paradijs 1989). 



6.1.2 Effects of radiation pressure on radial velocity 

curves 

The distortion of the secondary also affects the interpre- 
tation of radial velocity curves obtained from the phase- 
dependent Doppler shifts of spectral lines. In the presence 
of irradiation, the temperature of the heated side of the sec- 
ondary may be increased by a factor of 2 or more, depend- 
ing on the type of system and X-ray luminosity. This causes 
phase-dependent asymmetries in the shapes of individual 
spectral lines and shifts the effective velocity centre of the 
lines. If these effects are not properly taken into account, 
this leads to incorrect determinations of the radial veloc- 
ity amplitude (see Wade & Horne 1988; Shahbaz & Wood; 
Phillips, Shahbaz & Podsiadlowski 1999). The deformation 
of the Roche lobe introduces a further higher-order error. In 
order to illustrate the magnitude of these effects, we consider 
flux- weighted residual radial velocity curves calculated using 
the same system parameters as above (with Mi = 1.4 M©, 
Ma = IM0, Porb = 18.9 hr, = 2 x lO^^ergs"^). 

The lower panel of Figure 6 shows the residual radial ve- 
locity curves calculated with these parameters and an incli- 
nation of 40° for a model that includes only the variation of 
the surface temperature (including irradiation) and a model 
that includes the effects due to the temperature variation 
and the surface distortion. Here, the residual velocity was 




Figure 6. Top panel: The radial velocity curve of the centre of 

mass as a function of phase. Bottom panel: The residual radial 
velocity curves as a function of phase for the standard Roche lobe 
and radiation pressure model. The relevant parameters are given 
in the text. 



calculated according to 

J2r,e,^ Vr.d{r, e, 0, i)AF{r, 6, ct>, i) 



Kesid(4>, 



-VS^{^,i), 



AF(r,e,<A,$,i) 



(46) 



where Vrad{r,0,<j)j^j'i) is the radial velocity for a given or- 
bital phase angle $ and inclination i at a grid point (r, 0, 0), 
A/(r, 6, <j), ^, i) is the observed flux for a given surface ele- 
ment, and V^^{^,i) is the centre-of-mass radial velocity of 
the secondary. (The summation is only performed over the 
visible area of the star.) 

At phase 0, the observer sees the backside of the sec- 
ondary, and the radial velocity (and residual radial velocity) 
is zero. The radial velocity increases and becomes positive 
as the secondary moves away from the observer. However, 
the heated face of the secondary is now visible, and this 
shifts the 'effective centre' of the velocity towards the cen- 
tre of mass of the binary orbit. Hence the radial velocity 
is reduced and the residual radial velocity is negative. It 
falls to zero once more at phase 0.5, when the heated face 
of the secondary is directly in the observer's line-of-sight. 
After phase 0.5, the secondary moves back towards the ob- 
server and so the reverse effect is seen. The amplitude of 
these effects is substantially smaller (by about 30% or ~ 4 
km s^^) for the pressure modified surface compared with 
the standard Roche lobe. Again, this is a consequence of the 
decrease in the projected heated area perpendicular to the 
observer's line of sight, due to the deformation of the sur- 
face. This demonstrates that the distortion of the surface 
reduces the K correction that needs to be applied to ob- 
tain the true underlying radial velocity curve. This will be 
particularly important for systems like Sco X-1, where the 
measured X-ray luminosity (~ 1.8 x 10^* ergs^^; Kallman 
et al. 1998) indicates strong heating at phase 0.5. 
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7 CONCLUSIONS 

We have shown that irradiation of the secondary in a close 
binary can have a significant influence on the shape of irradi- 
ated stellar surfaces and on system configurations. We have 
constructed approximate numerical solutions of the equipo- 
tential surfaces of secondaries where external irradiation ef- 
fects are important to demonstrate the main physical effects. 
These generally scale with the luminosity of the radiation 
source and become dramatic whenever the ratio of radiation 
to gravitational forces approaches the critical value. The in- 
clusion of the radiation pressure term in the total effective 
potential leads to considerable deviations from the standard 
Roche geometry. Depending on the luminosity of the radi- 
ation source and the relative efficiency of absorption in the 
irradiated photosphere, the most important effects are: 

(1) The geometry of the equipotential surface in the 
vicinity of the irradiated region is altered, in particular for 
high X-ray luminosities. 

(2) For systems without accretion discs, the position of 
the irradiated inner Lagrangian (Li) point is shifted. This 
leads to an increase in the size of the critical Roche lobe. 

(3) When the irradiating luminosity increases above a 
critical value, the (irradiated) inner Lagrangian point ceases 
to be the Lagrangian point with the lowest potential. Then 
the critical potential will be set by the L2 point (or the L3 
point in HMXB's) and an outer critical configuration will 
be obtained. 

(4) In such outer critical configurations, mass loss oc- 
curs through the outer Lagrangian point. This may influence 
the evolutionary processes of systems such as binary pulsars 
or stars orbiting supermassive black holes in AGN. 

(5) In the case of mass loss purely from the L2 point, 
the associated change in angular momentum can lead to 
self-sustaining mass loss. This may be relevant for binary 
systems which show a steady decay in orbital period, and 
may lengthen the stable lifetimes of some HMXBs. Centau- 
rus X-3 would appear to be a possible candidate. 

(6) Due to the deformation of the stellar surface, 
lightcurves and radial velocity curves in LMXBs and 
HMXBs are modified. These effects need to be included for 
a reliable determination of system parameters. 

Finally, we note that while it has yet to be demon- 
strated conclusively that radiation pressure can effect the 
global structure of a star to the extent proposed in this pa- 
per, it is clear that the pressure force due to external ir- 
radiation can be of a considerable magnitude. No analysis 
of high-luminosity X-ray binary systems can be complete 
without an adequate treatment of these effects. 



APPENDIX 
CODE 



THE OPTICAL LIGHTCURVE 



Surface Geometry: 
The secondary is modelled using a surface modified by 
irradiation pressure with the assumptions of synchronous 
rotation and a circular orbit. The surface is then divided 
into grid points spaced equally in cos 6 and (j>. The accretion 
disc is modelled as described in S4.2.2. 



Von Zcipel's (1924b) theorem provides a relationship be- 
tween the local potential gradient and the local emergent 
flux / in a tidally or rotationally distorted star in radiative 
equilibrium, /rad = crT^ oc | Vfij. Consequently, the temper- 
ature at any point on the star is given by 

(1) 



pole 



\Vfl\{x,y,z) 



pole 



where Tpoie and |VO|poie are the temperature and potential 
gradient at the pole of the star. The gravity darkening 
exponent (5 has two values: 0.25 for stars with radiative 
atmospheres (von Zeipel 1924b), and 0.08 for stars with 
fully convective envelopes (Lucy 1967). 

Limb-Darkening: 

For each grid point on the secondary surface, the tempera- 
ture obtained from gravitational darkening is used to com- 
pute the monochromatic intensity according to Planck's re- 
lation: 



7(A,T) cx [exp (fec/fcAT) - 1]" 



(2) 



This intensity is modified using a standard linear limb- 
darkening law, 



I(fi,) = /(!){! — u + MC0S7}, 



(3) 



where 7 is the angle of foreshortening, and /i = cos 7. /(/n) 
represents the distribution of the emergent intensity which 
varies with the angle of foreshortening. The limb-darkening 
coefficients u are taken from standard tables computed 
from model atmospheres (e.g. Al Naimiy 1978). See Kopal 
& Kitamura (1968) and Avni (1978) for more discussions 
on these approximations. 

X-ray Heating: 

The X-ray flux that each surface element intercepts is given 
by 



fx-va^y{x,y,z) 



■ cose{x, y,z), 



(4) 



where Lx is the luminosity of the source, d the distance from 
the element to the source, and e the angle of foreshortening. 
Assuming total tliermalisation of the incident radiation, we 
can determine the modified temperature of each surface el- 
ement visible from the source (i.e. for which cos e > 0) from 



Tx-r!,y{x,y,z) 



'- pole 



\Vn\{x,y,z) 



4/3 



-I- 



|Vn|pole 
{l-W)U-r.yix,y,z) 



(5) 



where W is the albedo, and a the Stefan constant. 



We now have a grid of temperature values over the sec- 
ondary surface from which we can calculate the total lumi- 
nosity and flux. This calculation is repeated in steps of the 
phase angle. Optical lightcurves and flux-weighted radial ve- 
locity curves may then be calculated. 



Gravity-Darkening: 
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